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FOREWORD 


For  planning  purposes  and  questions  concerning  sensor  coverage, 
it  is  useful  to  know  the  analytical  ground  track  of  an  earth  satellite. 
This  report  provides  methods  for  determining  this  ground  track. 

This  technical  report  has  been  reviewed  and  is  approved. 

THOMAS  0.  WEAR,  Colonel,  USAF 

Director,  Space  Defense  Systems  Program  Office 

Deputy  for  Surveillance  and  Control  Systems 
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ABSTRACT 


This  report  derives  the  analytical  expressi 
satellite  ground  track  for  a  general  elliptical, 
circular-synchronous  orbit  under  the  assumptions 
conditions  and  in  the  absence  of  perturbations. 


ons  for  an  earth 
a  circular  and  a 
of  two-body 


SECTION  I 


INTRODUCTION 


The  purpose  of  this  paper  is  to  derive  parametric  formulas  for 
the  ground  track  of  an  artificial  earth  satellite.  Certain  assumptions 
are  made,  namely,  that  the  earth  is  a  sphere  and  that  two-body  orbital 
mechanics  hold,  i.e.,  the  satellite  undergoes  no  perturbations.  After 
deriving  the  general  expression,  two  particular  cases  are  considered. 
These  are  the  case  of  a  circular  orbit  and  the  case  of  a  circular 
orbit  of  24-hour  period  which  is  inclined  to  the  earth's  equator.  In 
these  cases,  time  is  eliminated  from  the  expressions  and  the  longitude 
of  the  sub-satellite  point  is  expressed  as  a  function  of  latitude. 
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SECTION  II 


GENERAL  CASE 

Consider  a  right  handed  inertial  x-y-z  coordinate  system  with 
origin  at  the  center  of  the  earth,  with  the  x-y  plane  coincident  with 
the  plane  of  the  earth's  equator  and  with  the  z-axis  intersecting  the 
north  pole.  Finally,  let  the  x-axis  point  towards  the  vernal  equinox. 
Also  consider  a  rotating  coordinate  system,  x^-y^-z^,  with  origin  at 
the  center  of  the  earth,  the  x^-y^  plane  rotating  in  the  x-y  plane 

A  A  A 

and  with  the  z^-axis  coincident  with  the  z-axis.  Let  i,  j  and  k  be 

A  A  A 

orthogonal  unit  vectors  in  the  x-y-z  system  and  i^,  and  be 
orthogonal  unit  vectors  in  the  x^-y^-z^  system.  Let  R  be  a  vector 
from  the  center  of  the  earth  to  the  satellite  and  let  the  satellite 
have  coordinates  (x,  y,  z)  in  the  x-y-z  system  and  (x^,  y^,  Zj)  in 
the  Xj-yj-Zj  system.  Finally,  let  the  period  of  the  rotation  of  the 
Xj-y^  plane  be  24  hours  so  that  its  angular  frequency 
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From  the  above  diagram. 


so  that 

^  A  A  A  ^  ^  A 

R  =  xi  +  yj  +  zk  =  x^  +  y]j]  +  z]k] 

A  _»  A  A  A  A  A  A 

x^  =  ij.R  =  xi^.i  +  yi  p j  +  zipk. 

Now,  the  angle  between  the  x  and  the  x^  axes  will  be  G^t  if  the  x^ 
axis  pierces  the  point  of  intersection  of  the  prime  meridian  and  the 
earth's  equator  and  t  is  the  time  in  mean  siderial  hours  from  when 
these  two  axes  coincide.  Therefore, 


A  A  A  A 

i  1 .  i  =  cosOJ^t  and  i  j. j  =  sinOJ^t. 

Also, 

A  A  A  A 

i  pk  =  0  since  k  =  ky 

Thus, 

A  A  ^  A 

x1  =  xcosUJ^t  +  ysinCJEt,  =  x j  ^ .  i  +  yj  y j  = 

-  xsinGJ^t  +  ycosOd^t,  and  z^  =  z. 

From  any  book  on  celestial  mechanics. 


where  a 

x  =  aP^  (cosE-e)  +  bQ^sinE 

y  =  aP  (cosE-e)  +  bQ.  sinE 

7  y  y 

z  =  aPz  (cosE-e)  +  bQzsinE 

is  the  semi -major  axis  of  the  satellite's  orbit,  e  is  the 

orbit's  eccentricity  and  E  is  the  eccentric  anamoly  of  the  satellite 
at  time  t.  E  is  given  by  Kepler's  equation 

E  -  esinE  =  2 TT  (t  -  T  ) 

—  ° 

o 
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where  Tq  is  the  time  of  perigee  passage  of  the  satellite  and  Pq  is 
the  orbital  period  of  the  satellite.  All  times  are  measured  in  the 
same  units  as  t.  The  semi -minor  axis  of  the  satellite  orbit 


b  =  a  (1  -  e2)2. 


P  =  cos  XL  cosco -  si  n  XI  sinGJcosi 
x 

Q  =  -cos-0-  sinOJ-  sin -O.  cos  cJcosi 
x 

P  =  sin-ft  cos +  cos-H-si  no^cosi 

y 

Q.  =~si n  -0-  si  n OJ+  cos  Si  cosOJcosi 


P^  =  s  i  n  s  i  n  i 

Q.  =  cos^>sini 
2 


where  .fLis  the  right  ascension  of  the  ascending  node  of  the  orbit 
relative  to  the  vernal  equinox, GJis  the  argument  of  perigee  of  the 
orbit  and  i  is  the  inclination  of  the  orbit  to  the  earth's  equator. 
Subs ti tuti ng. 


(cos-fl  sinGo  +  sinJl  cosGJcosi)  sin  EJ  cos  u>Et  + 


=  J^a  (sin/I  cosW+  cos  Si  si nCOcosi )  (cosE  -  e)  +  b  * 


2^  =  a  (sinwsini)  (cosE  -  e)  +  b  (coscjsini)  sinE. 
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Let  A  be  the  longitude  of  a  point  on  the  earth's  surface  measured 

from  the  prime  meridean  eastward  from  0°  to  360°.  Let/^  be  the  latitude 

of  the  same  point  measured  north  of  the  equator  from  0°  to  +90°  and 

o  •  o 

south  of  the  equator  from  0  to  -90  .  If  this  point  is  the  sub-satellite 
point  of  the  ground-trace,  then  the  satellite's  location  on  the  celestial 
sphere  will  have  the  same  latitude  and  longitude  as  the  sub-satellite 
point,  since  the  celestial  sphere  is  concentric  with  the  sphere  of  the 
earth.  Also,  let 

R  =  1R1. 


Then, 


and 


and 


=  Rcos  A  cos y3  ,  y^  =  Rsin  A  cos^  , 

Zj  =  Rsin^  so  that  cos  A  cosfl  =  ^  (xcosCJ^t  + 

R 

ysi  nCJEt), 

sin  A  cos/J  =  J_  (ycosOJpt  -  xsintJpt), 

'  R 


sin/?  =  z  . 

R  * 

substi tuti ng 

R  =  a  (1  -  ecosE) 


and  using 

b  =  a  (1  -  e2)K 

tan  A  =  £  [ ( si n il.  cosW+  cos./!  sincjcosi)  (cosE  -  e)  + 

2  1  , 

(1  -  e  )2  (cosI\  cosOJcosi  -  sin.fl  sinu>)  sinEjx 
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coscJEt  -  [(cos  XI  COSOJ-  sinil  sinGJcosi )  (cosE  -  E)  - 
(1  -  e^)2  (cosJl  sinCL>+  sin  SX  cos  GOcos  i )  sinE^J  sinGJ^t^}* 
[  [(cos /I  cost*)-  sin SL  sinujcosi )  (cosE  -  e)  -  (1  -  e^)"1-  >. 

(cosil  sinu>+  si  nil  cosodcosi )  sinE*J  cos^t  + 

r  2  — 

^(si nil  cosuj+  cosil  sintucosi)  (cosE  -  e)  -  (1  -  e  )2  * 

(sinjfL  sinw-  cos  J1  coscocosi )  si  nE^]  si  n  GJ^t}  ~  ^ . 

r  2  — 

sin/?  =  [(sinojsini)  (cosE  -  e)  +  (1  -  e  )2  (cosojsi ni  )* 


s  i  nE  ^  [ 


1  -  ecosE 


The  quadrants  of  /(  and^  are  easily  determined  from  the  following 
considerations.  must  be  between  +90°  and  -90°.  If  sin^*  is  +,  the 
ground  track  is  in  the  northern  hemisphere,  andy^is  between  0°  and  90°. 

If  sin^?  is  the  ground  track  is  in  the  southern  hemisphere  and^^is 
between  -90°  and  0°.  The  quadrant  for  X  is  determined  in  the  conventional 
way  be  considering  the  signs  of  and  y^;  i.e., 
tan X  =  y,  • 

T\ 
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SECTION  III 


CIRCULAR  ORBIT 


Two  special  cases,  often  useful  for  various  planning  purposes, 
are  discussed  next.  The  first  of  these  is  the  circular  orbit.  Here, 
e  =  o.  Also,  E  =  (x>st  where  c«>s  is  the  angular  frequency  of  the 
satel 1 i te  and 

( jJ .  =  2 £_. 

P 

o 

Since  the  argument  of  perigee  is  arbitrary  for  a  circular  orbit, 
we  may  set  o,  i.e.,  the  line  of  apsides  coincides  with  the  line  of 
nodes.  Also,  R  =  a  for  all  t.  Substituting  these  values, 

tan  A  =  £  (si  n  A  cos  cjgt  +  cos  fl  cosisinodgt)  costO^t  - 
(cosXL  cos6Jst  -  si  nil.  cosisincJ^t)  sinCJ^t^jj^  (cosiT  * 
cosCJgt  -  sin  X"LcosisinCJst)  cosOOgt  +  (si  n XL  cos  (aM  + 
cos  XI  cosisin  G.'st)  sin6J^t"j  '  j 
sir^=  si  nisincJgt. 

Since  from  this  relation  for  sin^y  it  follows  that 

.  -1 


t  =  J _ 

OJ 


sin'  /  sin^  , 
\  sim  / 

we  can  write  ^explicitly  as  a  function  of^  : 

i/V=  £^sinjQ  cos  f  si, 


tan , 


cos  i 


h 


cjc 
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si n SI  ctni sin 


sin''  (f^) }  * 

-1 


It  should  be  noted  that  for  a  giveny?  there  will  be  two  values 
of  A  per  or bi t, 
one  for 


cosJT.  ctni sin^\sin(-Gk>£  sin"' 

/  te 


si  rti 


and  one  for 


t  =  J T  s  i  n  '  f  s  i  n  /S  \  +  Tf  j  , 

oj$  L  \  sirfi  J  2  J 

the  sign  of  +  7 T  corresponding  to  the  sign  of  sin  ' 

2 


A  further  simplification  will  result  if  the  x-axis  is  rotated  such 
that  it  coincides  with  the  ascending  node,  i.e.,  such  thatX^.=  0.  Then 
the  longitude  of  the  ascending  node  relative  to  the  prime  meridian,  XI  Q> 
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must  be  used  in  constructing  plots.  This  simplification  is  utilized 
in  the  next  section. 
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SECTION  IV 


CIRCULAR  SYNCHRONOUS  ORBIT 

The  next  interesting  special  case  occurs  when  the  orbit  is 
circular  with  the  period  the  same  as  that  of  the  earth,  i.e., 

H  =WE> 

and  when  the  direction  of  motion  is  the  same  as  the  direction  of  the 
earth's  rotation.  Here  we  may,  for  convenience,  let  the  x-axis  pierce 
the  ascending  nodal  point  of  the  orbital  plane,  that  is,  setfl=  0. 

We  again  get  ^  as  a  function  of^  by  the  following  method: 


-  sinCJ^t  cosC*^t  = 

(cosi  -  1)  sinCJ^t  coscJ^t; 
sinCJ^t  =  sin/*?  ; 


sTni 


therefore,  sin./\  cos^  = 


sin 2/S  )2 


By  knowing  the  inclination  of  the  orbit  of  such  a  satellite,  we 
may  easily  compute  its  ground-trace  and  from  this  deduce  the  satellite’s 
earth  coverage  at  any  instant  of  time. 

It  is  necessary  to  reiterate  that  this  ground-trace  will  not  be 
relative  to  the  vernal  equinox  but  will  be  relative  to  the  ascending 
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node  of  the 
stati onary. 


orbit.  Also,  when  i  =  0,  the  orbit  is  synchronous 
i .  e . ,  y?  =  \  =  0  for  all  t . 
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SECTION  X 


CONCLUSIONS 

In  many  instances,  where  only  a  rough  approximation  is  required, 
the  two-body  ground-trace  will  suffice.  The  analytical  formulas 
presented  will  rapidly  provide  the  trace  and  a  particularly  interesting 
case,  the  circular  synchronous  orbit,  is  especially  amenable  to  hand 
calculation. 
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